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We prove an upper bound for the convergence rate of the homogenization limit e — > for 
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justify the formal homogenization asymptotics obtained by us earlier by proving an upper 
bound for the convergence rate (a corrector estimate). The main ingredients of the proof 
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understanding of the analysis of the limit two-scale problem. 
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1 Introduction 



We study the averaging of a system of reaction-diffusion equations with linear trans- 
mission condition posed in a class of highly heterogeneous media including areas of low 
and high diffusivity. Our aim is twofold: on one hand, we wish to justify rigorously the 
formal asymptotics expansions performed in 30 , while on the other hand we wish to 



understand the error caused by replacing a heterogeneous solution by an approximation 
(averaged) estimate. To this end, we prove an upper bound for the convergence rate of 
the limit procedure e — > 0, where e is a suitable scale parameter (see section 2.1 for 
the definition of e). The materials science scenario we have in view is motivated by a 
very practical problem: the sulfate corrosion of concrete. We refer the reader to [3] for a 
nice and detailed description of the physico-chemical scenario and to 17 for the formal 



(two-scale) averaging of locally-periodic distributions of unsaturated pores attacked by 



sulphuric acid. The reference 18 contains the rigorous proof of the limiting procedure 
e — >• treating the uniformly periodic case of the same corrosion scenario. In p], the 
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main working tools involve the two-scale convergence concept in the sense of Nguetseng 
and AllaircQ combined with a periodic unfolding of the oscillatory boundary. Here, we 
use an energy-type method. 

The framework that we tackle in this context is essentially a deterministic one. We 
assume that a distribution of the locally-periodic array of microstructure^ is known a 
priori. We refer the reader to [6] for related discussion of continuum mechanical descrip- 
tions of balance laws in continua with microstructure. At the mathematical level, we 
succeed to combine successfully the philosophy of getting correctors as explained in the 
analysis by Chechkin and Piatnitski [7] with the intimate two-scale structure of our sys- 
tem; see also [8j [9] for related settings where similar averaging strategies are used. For 
methodological hints on how to get corrector estimates, we refer the reader to the analy- 
sis shown in 14 for the case of a reaction-diffusion phase field-like system posed in fixed 
domains. For a higher level of discussion, see the monograph 10 which is a nicely writ- 
ten introductory text on homogenization methods and applications. It is worth noting 
that the periodic homogenization of linear transmission problems is a well-understood 
subject (cf. [l] [l5j [l9j [28] , e.g.), however much less is known if one steps away from the 
periodic setting even if one stays within the deterministic case. If the microstructures are 
distributed in a suitable random fashion, then concepts like random fields (see [4j or [5] 
and references cited therein; 21 and intimate connection with the stochastic geometry 
of the perforations) can turn to be helpful to getting averaged equations (eventually also 
capturing new memory terms), but corrector estimates seem to be very hard to get; cf., 
e.g., [2j. Disordered media for which the stationarity and/or ergodicity assumptions on 
the random measure do not hold or situations where at a given time scale length scales 
are non-separated are typical situations that cannot be averaged with existing techniques 
(so it makes no sense to search here for correctors). 

This paper is organized as follows: In Section[2]we introduce the necessary notation, the 
microscale and the two-scale limit model. Section [3] contains the main result of our paper 
- Theorem 3.1 In Section [4] we introduce the technical assumptions needed to obtain the 
correctors. At this point, we also collect the known well-posedness and regularity results 
for both the microscopic model and the two-scale limit model. The rest of the paper 
consists of the proof of the convergence rate and is the subject of Section [5| 



2 Statement of the problem 

In this section we introduce our notation, the microscale model and the two-scale limit 
problem. We start off with the definition of the locally periodic heterogeneous medium. 
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1 For an introduction to two-scale convergence, see 

2 Here, we deviate from the purely periodic setting. It is worth noting also that the assump- 
tion of statistically homogeneity of distribution of microstructures, which is a crucial restriction 
of stochastic homogenization, does not cover all possible configurations of locally-periodic ge- 
ometries. Also, as we formulate the working framework, we cannot treat neither randomly placed 
microstructures nor stochastic distributions of microgeometries. Consequently, the exact con- 
nection between these two averaging techniques is hard to make precise. 
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2.1 Definition of the locally-periodic micro medium 
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Figure 1. Schematic representation of a locally-periodic heterogeneous medium. The cen- 
ters of the gray circles are on a grid with width e. These circles represent the areas of 
low diffusivity and their radii may vary. 

We consider a heterogenous medium consisting of areas of high and low diffusivity. 
The medium is in the present paper represented by a two dimensional domain. For the 
definition of the locally-periodic medium with inclusions, we follow in main lines (7). We 
denote the two dimensional bounded domain by fl C K 2 , with boundary T. Denote 

J £ :={jeZ 2 |dist(ej,r)>e\/2}, 

U := {y E K 2 | - 1/2 < Vi < 1/2 for i = 1,2}). 

A convenient way to parameterize the interface r e between the high and low diffusivity 
areas, is to use a level set function, which we denote by S e (x): 

ier £ » s e (x) = o, 

Since we allow the size and shape of the perforations to vary with the macroscopic variable 
x, we use the following characterization of S e : 

S e (x) := S(x,x/e), (2.1) 

where S : fl x U — > R is 1-periodic in its second variable, and where S is independent of 
e. We assume that S(x, 0) < const. < and S(x,y)\ y& gu > const. > for all x G fi, so 
that the areas of low diffusivity in each unit cell do not touch each other. We set 

Q) :={xee(U+j)\S(x,x/e) < 0}, 
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and introduce the area of low diffusivity Oj as follows: 

m U Qp 

and the area of high diffusivity fl e h as: 

Q e h ■.= n\nf. 

We call a medium of which the geometry is specified with a level set function of the type 
that is given in (2.1 1 a locally periodic medium |7j. In Fig. [I] a schematic picture is given 



of how such a medium might look like for a given e > 0. Note that by construction the 
area of high diffusivity tt% is connected and the area of low diffusivity fif is disconnected. 
The interface between high and low diffusivity areas T e is now given by r e = d£lf and 
the boundary of Vl\ is given by <9S1| = T c U T. 

Furthermore, we use the notation Q,\ := f2\ (J - e Jr _ (e(U +j)), and we introduce for later 
use the smooth cut-off function Xe(%) that satisfies < Xe{ x ) < 1; Xe( x ) = if x € Q\ 
and Xe(x) = 1 if dist(x,fif) > dist(r e ,Of). Moreover, e|Vx £ | < C and e 2 |Ax e | < C, with 
C independent of e, and also 

\\l-Xe\\mn) <e 1/2 C, 

l|Vx e ||L»(n) < e" 1/2 C, (2.2) 

||Ax e ||i»(n) <e- 3/2 C, 

with C again independent of e (see e.g. [7], [14] ) . 

In addition, we need to expand the normal v e to T e in a power series in e. This can be 
done in terms of the level set function S t , which we assume to be sufficiently regular so 
that all the following computations make sense (see assumption (Bl) in section [4j: 

VS'Qe) _ VS(x,x/e) V x S+lV y S 

\VS^x)\ \VS(x,x/e)\ |V x 5+iV„S| " [ ' 

First we expand |ViS e |. Using the Taylor series of the square-root function, we obtain 

|V5 e | = ±\V y S\ + 0(e°). (2.4) 



In the same fashion, we get 



where 



t/ = v Q + evi + 0(e 2 ), 



\VyS\ 



and 



(y x s-v y s) v y s 

\VyS\ \VySf \VySl ^ 

If we introduce the normalized tangential vector To, with To J- vq, we can rewrite v\ as 

^ = T0 ^sT- (2 - 6) 
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Later on we will also use the notation 

v\ =£-V^o) =2^1 + 0(e). 
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2.2 Microscopic equations and their upscaled form 

We focus on the following microscopic model 

[d t u e = V • (D h Vu e - q e u £ ) mn e h , (2.7) 

{d t v e = e 2 V • (DiVv e ) in Of, (2.8) 

= eV • (DiVv e ) , , , 

V ; onr £ , (2.9) 

U e (x,t) =u b (x,t) on T, (2.10) 
I u t (x,0) = u[{x) in 




i> e (a;, 0) = v\ (x) in £1 



(2.11) 



Here we denoted the tracer concentration in the high diffusivity area by u e , the con- 
centration in the low diffusivity area by v e , and the velocity of the fluid phase by q e . 
Furthermore, denotes the diffusion coefficient in the high diffusivity region, D; the 
diffusion coefficient in the low diffusivity regions, v e denotes the unit normal to the 
boundary T e (t), where u& denotes the Dirichlet boundary data for the concentration 
u e , and where u\ and vl are the initial values for the active concentrations u e and v e , 
respectively. 

In order to formulate the upscaled equations and obtain closed-formulae for the effec- 
tive transport coefficients, we use the notations 

B(x):={yeU:S(x,y)<0}, (2.12) 

and 

Y(x) :=U -B{x), (2.13) 
and we define the following cc-dependent cell problemj^] 

A y Mj(x, y) = for all x € fi, y € Y(x), 

u ■ V y Mj(x, y) = — uq ■ ej for all iGllje dB(x), (2.14) 
Mj(x,y) y-periodic, 

for j 6 {1,2}. To ensure the uniqueness of weak solutions to this cell problem, suitable 
conditions on the spatial averages of the cell functions need to be added. 

The solution to this cell problems allows us to write the results of the formal homog- 
enization procedure in the form of the following distributed-microstructure (two-scale) 



3 See, e.g. 
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, for the role of cell problems in formulating averaged equations. 
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model 

d t v (x, y, t) = DiA y v (x, y, t) for y E B(x), x E f2, 

d t {e(x)u + J lyl<r{x) v dy) = (2.15) 

V x ■ (D(x)V x u - qu Q ) for x E il, 



\v (x,y,t)=u (x,t) for y E dB(x), 
I Uo(x,t) = Ub(x,t) for ier, 

{uq(x,0) = uj(x) for x E 0, 

v (x,y,Q) = vi[x,y) for y E B(x), x E il. 

where the porosity 9{x) of the medium is given by 

6{x) := \Y(x)\=l-\B(x)\, 

and where the effective diffusivity tensor T>(x) is defined by 



with M = (Mi,M 2 ). 



(2.16) 
(2.17) 



2>(a:) :=D h / (7 + V„M(a:,v)) dy, 

JY(x) 



3 Main result 



Recall that (u e , v e ) is the solution vector for the micro problem and (uq, vq) is the solution 
vector for the two-scale problem. We introduce now the macroscopic reconstructions^] 
Uq , Vq , u\ , which are defined as follows 

Uq(x, t) := u n (x, t) for all x E (3-1) 
Vq(x, t) := Vq(x, x/e, t) for all x E Slf, (3.2) 
u\(x,t) := ul{x,t) + eM(x,x/e)Vu%{x,t) for all x E Q. e h , (3.3) 

In the same spirit, we introduce the reconstructed flow velocity q^x) = q(x) for all 
x E Sl e h and the corresponding reconstructions u e 0I ,Vqj for the macroscopic initial data 
ui and vi, respectively. 

The main result of our paper is stated in the next Theorem. The applicability of this 
result confines to our working assumptions (Al), (A2), (Bl), (B2), and (B3) that we 
introduce in Section [4] 

Theorem 3.1 Assume (Al), (A2), (Bl), (B2), and (B3). Then the following conver- 
gence rate holds 

\\ u t - Uq\\l^{I,L 2 (Q%)) + \\Ve - WolU°°(/,i 2 (on) + 

IK - mU~(/,HM^)) + e IK - Vq\\ l °° (i, < c Vt, (3-4) 
where I = (0, T], and where the constant c is independent of e. 



We borrowed this terminology from 14 . Note however that the concept of reconstruction 



operators appears in various other frameworks like for the heterogeneous multiscale method 13 
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The remainder of the paper is concerned with the proof of Theorem 



3.1 
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4 Technical preliminaries 
4.1 Function Spaces. Assumptions. Known results 

4.1.1 Functional setting 

As space of test functions for the microscopic problem, we take the Sobolev space 

H l {n e ;d£l) := {tp G ff x (O e ) with <p = at dtt} . 

Since the formulation of the upscaled problem involves two distinct spatial variables 
x € and y £ B{x) (with B(x) C fi), we need to introduce the following spaces: 

Vx:=H^(n), (4.1) 

V-2 ■= L 2 (ft; H 2 (B(x))), (4.2) 

ffi:=i§(fi), (4.3) 

iT 2 :=L 2 (r»;L 2 (B(x))). (4.4) 



The spaces H 2 and V2 make sense, for instance, as indicated in 26 



4.1.2 Assumptions for the microscopic model 

Assumption (Al): We assume the following restrictions on data and parameters: Take 
D h ,Di e]0,oo[, u{ £ H x {Q? h ), v{ e if 1 (Of), and w fc e i/ 1 (7; if 3 (r)). Furthermore, we 
assume 

lk-«ojlk<-(nj) = 0(e*) with ^ > l - 
H ~ «o/IU~(n f ) = <V 2 ) with 2 > i 

Assumption (A2): q t E H 2 ^; R d ) n L°°(0^ M d ) with V • <? £ = a.e. in £2£. Addition- 
ally, we assume that 

Ike - IqWl'W*') = 0{e^) with 9 3 > i (4.5) 

Remark 4.1 If one wishes to replace q e with the stationary Stokes or Navier-Stokes 
equations, then a few additional things have to be taken into account. One of the most 



striking facts is that the exponent 9 arising in (4.5) seems to be restricted to #3 = 
|. Consequently, this worsens essentially the convergence rate; see, for instance 



24 



(Theorem 1) for a discussion of homogenization of the periodic case. 



4.1.3 Assumptions for the two-scale model 

Assumption (Bl): The level set function S : f2 x U — > K is 1-periodic in its second 
variable and is in C 2 ( fl x IX). 
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Assumption (B2): We assume the following restrictions on data and parameters: 

9, d e £+ (0) nH 2 (Q), 

q £ H 2 (n;R d ) n L°°(ri;R d ) with V • q = a.e. in O, 

u b e Lf(n x i) n h 1 ^; # 3 (r)), 

d t Ub < a.e. (x, i) G Cl x I, 

u/ e L^°(0)niJinif 2 (0), 

^/(x, •) e n ff 2 for a.e. lefl. 



Assumption (B3): 



^ :=H x (n). 

H 2 :=£ 2 (0;tf 2 (£?(ar))), 

Vi :=# 9 2 (ft), 

F 2 :=L 2 (Sl;H 3 (B(x))). 



(4.6) 
(4.7) 
(4.8) 
(4.9) 



Remark 4.2 Following the lines of 26 and 29 , Assumption (Bl) implies in particular 



that the measures |9B(x)| and |B(x)| are bounded away from zero (uniformly in x). 
Consequently, the following direct Hilbert integrals (cf. [12] (part II, chapter 2), e.g.) 



U ( X )\\h(dB(x)) < °°} 



L 2 {VL- H l {B{x))) := {u e L 2 (ft; L 2 {B{x))) : V y u e L 2 (ft; L 2 {B{x)))} 
L 2 (Q; H l (dB(x))) := {u : fi x dB(x) -> K measurable such that 

are well-defined separable Hilbert spaces and, additionally, the distributed trace 
7 : L 2 {n-H 1 {B(x))) -> L 2 (fl, L 2 (dB(x))) 

given by 

7 u(x,s) := (7 x t/(x))(s), x G fi,s G dB(x), u € L 2 (Q; H 1 (B(x))) (4.10) 



is a bounded linear operator. For each fixed x £ J7, the map 7^, which is arising in (4.10), 



25 



is the standard trace operator from H 1 (B(x)) to L 2 {dB{x)). We refer the reader to 
for more details on the construction of these spaces and to [27] for the definitions of their 
duals as well as for a less regular condition (compared to (Bl)) allowing to define these 
spaces in the context of a certain class of anisotropic Sobolev spaces. 

For convenience, we also introduce the evolution triple (V, H,V*), where 

V:={(0,V) G Vi x V 2 \4>(x) = ^(x,y)forxe fl, y e dB(x)}, (4.11) 
U := R x x H 2 , (4.12) 



4.1.4 Analysis of microscopic equations 

Definition 4.3 Assume (Al), (A2) and (Bl). The pair (u e ,v £ ), with u e = U e + u b and 
where (U e ,v e ) € H 2 (fl^; dti) x H 2 {Vt\), is a weak solution of the problem (P e ) if the 
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transmission conditions (2.9) are fulfilled and the following identities hold 



d t (U e +u b )<j)dx+ / (D h V(U £ +u b ) - q e (U e +u b )) ■ Vfidx = 

-J v e • (e 2 Di\7v t )cj)ds, (4.13) 
/ d t v £ ip dydx + / e 2 D t Wv e ■ Vip dx = / v € ■ (e 2 Di\7v e )tp dsdx, (4.14) 
for all (0,-0) e H^n^dfl) x iP(Of) and t € I. 



Theorem 4.4 Assume (Al), (A2) and (Bl). Problem (P t ) admits a unique global-in- 
time weak solution in the sense of Definition \J~. 



Proof Since we deal here with a linear transmission problem, the proof of the Theorem 



can be done with standard techniques (see 16 , e.g.). □ 



(P){ 



4.1.5 Analysis of two-scale equations 

This section contains basic results concerning the well-posedness of the two-scale problem, 
which we reformulate here as: 

'9(x)d t u - V x ■ (D(x)V x u - qu ) = - SdB{x) ' (AV„«o) da in 0, 
d t v — DiAyVo = in B(x), 

uo(x, t) = vq(x, y, t) at (a?, y) G fi x dB(x), 

uo(x, t) = u b {x, t) at x € 90, 

Uo(x, 0) = Ui(x) in O, 

v (x,y,0) = vi(x,y) at (x,y) e x B(ie). 

Before starting to discuss the existence and uniqueness of weak solutions to problem 
(P), we denote U := u — u b and notice that U = at 90. 

Definition 4.5 Assume (Bl) and (B2). The pair (it, v), with u = U + u b where ([/, v) £ 
V, is a weak solution of the problem (P) if the following identities hold 



6d t (U + u b )(/)dx+ I (DV X (U + u b ) - q(U + u b )) -V x (j)dx = 

/ / v ■ (Di\7 y v)(j)dadx 

JQ JdB(x) 

(4.1 

dtvipdydx+ / / -D;V y • V y ip dydx — / / v§ • (DiV y v)(f>dadx, 

JnJB(x) JnJdB(x) 

for all (0, V>) € V and f e J 



nJB(x) JnJB(x) JnJdB(x) 

(4.16) 
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Proposition 4.6 (Uniqueness) Assume (Bl) and (B2). Problem (P) admits at most one 
weak solution in the sense of Definition 14.51 



Proof Since Problem (P) is linear, the uniqueness follows in the standard way and 
can be done directly in the x-dependent function spaces: One takes two different weak 
solutions to (P) satisfying the same initial data. Testing with their difference and using 
the Gronwall's inequality conclude the proof. □ 

Theorem 4.7 Assume (Bl) and (B2). Problem (P) admits at least a global-in-time weak 
solution in the sense of Definition \4-5\ 



Proof See the proof of Theorem 5.11 in 30 



□ 



To get the correctors estimates stated in Theorem 3.1 we need more two-scale regularity 
for the macroscopic reconstruction of the concentration field vq. We state this fact in the 
following result: 



Lemma 4.8 (Additional two-scale regularity) Assume (Bl) and (B2). Then 

v* eL 2 (I;H 2 (Sl;H 2 (B(x))))). (4.17) 



Proof The proof of this result is similar to the regularity lift proven in Claim 5.10 in 30 



The main ingredients are fixing of the boundary and testing with difference quotients in 
the weak formulation posed in fixed domains. We omit to show the details. □ 



Theorem 4.9 (On strong solutions) Assume (Bl), (B2), and (B3). Problem (P) admits 
one global-in-time strong solution. 



Proof Under the assumptions (Bl) and (B2), Theorem |4.7| guarantees the existence of 
global-in-time weak solutions. Relying on the assumption (B3), we can lift the regularity 



until getting a strong solution. A similar calculation is done in 16 (Theorem 5, pp. 
360-364). In particular, (B3) allows for a regularity lift such that 

l|V9 t uollL 2 (/xn) + l|V<9 t uo||L2( /x o xy) < c. 
For our purpose, we only need 

\\Vd t u £ \\ L 2 iIxnv <c, (4.18) 

IKIU s (i ; u-3(n«)) < c (4.19) 



where Uq is the macroscopic reconstruction of Uq. Quite probably (4.19 1 could be relaxed 
to ||wolli 2 (/;-H 2+8 (f2j )) f° r somc ^ > 0, but we don't address here this issue. Since we 
rather wish that the reader focusses on our strategy of getting the correctors, we omit to 
show the proof details for this regularity result. □ 
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In this section, we give the proof of the main result of our paper, i.e. of Theorem |3.1[ The 
proof uses the auxiliary results stated in the lemmas below. They mainly concern integral 
estimates for rapidly oscillating functions with prescribed average; Related estimates can 
be found, for instance, in [7] and Section 1.5 in [10] . 

Lemma 5.1 Assume the hypothesis of Theorem 3.1 to hold. Then 

/ f V* • ((/ + \7 y M)X7 x u ) - -^V, ■ / (/ + V y M)V x u dy) dy = 

JY(x) V °{ x ) JY(x) ' 



Proof We compute 



/ + V y M)V x u Q da. 

JdB(x) 



I ( V x ■ ((I + V y M)V / (/ + V y M)V x u Q dy) dy = 

Jy(x) v V\ x ) Jy(x) ' 

/ V a • ((/ + VyM^viio dy - V x • / (I + V y M)V x u dy. 

JY(x) JY(x) 

Reynolds's transport theorem (see for instance 20 ) gives 

V x - (I + V y M)V x u Q dy= I V x - ((I + V v M)V x u ) 

JY(x) JY(x) 



dy 



\VyS\ 



(I + W y M)V x u da. 



By the boundary condition in (2.14), we have 



vq(I + V v M) = 



VyS 

\VyS\ 



(I + VyM) = on dB{x), 



so that we can write, using (2.5), 
V S 

" (7 + V y M)V x u = 



W X S V x S-VySVyS 



\VyS\ 



(I + VyM)V 



JVyS\ \VyS\* \V y S\ f 

= !/!•(/+ V y M)W x UQ 

on dB(x). Combining the above expressions proves the conclusion of this lemma. 



□ 



Lemma 5.2 Assume the hypothesis of Theorem 3.1 to hold. Let Q(x, y) £ L 2 (fl; L 2 (B(x))) 
andp € L 2 (Q; L 2 (dB(x))). Furthermore, suppose that fyr x \ Q(x, y) dy — Jg Y r x ) P( x i v) 
Then the inequality 



/ 

isi: 



Q(x, x/e)(f>{x) dx — e / p(x, x/e)4>(x) ds 



< Cell 



lff l («£) 



holds for every <j> £ dil). The constant C does not depend on the choice of e. 
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Proof The problem 

A y ^(x,y) = Q(x,y) in Y(x) 
v ■ =p(x,y) on dB(x) 

has a 1-periodic in y solution that is unique up to an additive constant. We multiply the 
first equation above with (f> and integrate over D, e h and use the thus obtained equality to 
get: 



Q(x, x/e)(j)(x) dx — e / p(x, x/e)(j){x) ds 
A y 1 i'(x, y)\ y =*4>(x) dx — e / p{x 1 x/e)ds 
ej (v x [V y ^(x,y)\ y=li )-V x Vy9(x,y)\ y= ^(f)(x)dx-e I p(x,x/e) 
e {v Q + ev{)- V y ^{a,y)\ y= *(j)(x)ds- e I V y ^\ y= ^\7 x <j)(x) dx 

h 

-e VzVy^x, y)\ y =*(f)(x) dx - e / p(x,x/e)ds 



ds 



v\ ■ V y *(a, y)\ y= *(j)(x) ds 



+ € 



^x^ y ^{x,y)\ y= ^(t>(x) dx 



/ V y ^\ y= *V x (j)(x) dx 
<eCU\\ m{nV 



The lemma is now proved. 



□ 



The last auxiliary lemma is a special case of Lemma 4 in [7; , and therefore we will 
state it here without proof. 

Lemma 5.3 Let LI e be a subset of {x € fl | dist(x, dfl) < ^e}. Then the following 
inequality 



V x uq(/) dx 



< Ce^ 2 U\\ m{nv 



holds for all 4> £ H l {Sl € h ; dil). The constant C does not depend on e. 



Proof (of Theorem 3.1) We define 



z e {x, i) :— uq(x, t) + eXt(x)M{x, x/e)Vuo(x, t) — u e (x, t), 
w e (x, t) := vq(x 1 x/e, t) — u e (x, t). 

By, e.g., Theorem 4 in [23], the functions M(x,x/e) and v (x, x/e, t) are well-defined 
functions in and L 2 (I \ H 1 (Vlf)) , respectively, and furthermore, we know by con- 



Corrector estimates of a locally-periodic scenario 
struction that there exists a C > such that 



13 



IW|z,2(/;ifi(n«)) < C, (5.1) 
IKI|z,2(/;Hi(n«)) < eC 1 , (5.2) 

where the constant C is independent of the choice of e. In the following we will use the 
notation u\(x,y,t) = M(x,y)VuQ{x,t). We compute: 

Az e (x,t) =Au {x,t) + e Xe A x Ui(x, y, t)\ y= z + 2\ e S/ x ■ V y u\{x, y, t)\ y= ^ 
+ ^Xe\ui{x,y,t)\ y= <i + eA Xe Ui(x, y, i)| y =f - Au e (x,t) 
+ 2eVx e • V x u 1 (x,y,t)\ y= <i + 2\7 Xe ■ V y ui(a;, y, i)|y=f 
AxW e (i, =A x u (a;,i) + e~ 1 V a; • V^uo^, 2/, i)|y=f 

+ e _1 Vj, • V x u (a;,y,i)|j /= | + -^A y v (x, y, t)\ y= ^ - A x w e (a;, y, t)\ y= ^ . 

We use that 

9(x)d t u - V x • (2?(a;)V a; Mo - g"o) = / • (D[V y v ) dain fl 

JdB(x) 



and A y ui(x, y) — and <9 t u £ = V • (DhVu e — <7 e w £ ) to obtain 

D h A x z e - <9 t z £ =D h A x u + eXeD h A x ux + 2\eD h V x ■ V y ui + eD h AxeU\ 
+ 2eD h W X e ■ V^ii + 2D h W X e ■ V y ui 

- \ fix ■ (D(x)V x u - qu ) - I v - (DiW y v ) do) - t X Au x 

V V JdB(x) ' 

- V • (q t u £ ), 

e 2 DiA x w iL - d t w e =e 2 DiA x v + eDi{V x ■ V y v + V v ■ W x v ). 
On r c we have 



i/ e ■ Vz £ — — v e ■ Vu e + v e ■ V ' x u + ev e ■ V ' x ui + u e ■ V y u\ 

= - e 2 ^-v e ■ Vv £ + f £ • V x u + eu e ■ V ' x ui + v e ■ 
L>h 

v e ■ Vw £ — — v e ■ Vti £ + v e ■ V x v + e _1 i/ e • V y v a . 
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Now, we multiply with <f> and integrate by parts to get 



/ d t z e (f)dx+ I (-D/jV z t ■ \7(j) + q e ■ V z f jfi) dx = e / xA u i4>dx — / D^^UQcj) dx 
- e / DhX^xU^dx - 2 / B h x^S x ■ VyU^dx 

+ / \(v x -{T>{x)V x u Q )- ( v -{D l V y Vo)da)<t>dx (5.3) 
-D;z/ • Vv e (pds — Dh \ v e ■ V x uq4> ds 



eDh J v e ■ V x ui(f> ds — Dh J v e ■ S7 y ui(f> ds 
/ {-q- q t ) -Vuo4>dx - e j XeUiq e ■ Vfydx, 



(5.4) 



and 



/ dtw e Tp dx + e 2 / DiX?w e ■ Vt/> dx — — e 2 / DiA x voipdx 
Juf Jnf Jn c t 

— e Di(S7 x ■ \7 y vo + V y ■ \7 x vo)ip dx ~ e 2 / Div e ■ Vv e ip ds 

+ e 2 J D t v e ■ \7 x v ^ ds + e J D t v e ■ \7 y v ^ ds. (5.5) 

We take into account the identity 

(V„ • V !E « 1 (a;,y))| 1 , =a ./ e = tV x ■ {V x ui{x, x/e)) - e(A x ui(x,y))\ y=x/e , 
which gives 

eD h / v e ■ V x u 1 \ y=x / e z e ds —eD h / V x Ui \ y - x / € ■ V(x £ ■ xUl\y= x /e) dx 

JF' Jni 

h 

= eD h / V x Mi • V(XeZe) + XeZe^xUl dx 

Jat 



+ D h XeZe^y • V^Ui cfc, 

and also the boundary condition Vq ■ V y Ui = —v Q ■ V x u Q for y 6 dB{x), which holds 



also on r e , we add the two equations (5.3) and (5.5), substitute <p = z f _ and ip — w e and 



obtain 
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ld t \\z e \\ 2 L2{nv + ^KHi^n.) + f (D h Vz e • Vz e + q e ■ Vz € ) dx + e 2 D l \\Vw e \\ 2 L2{nD = 

Jn l, 

-[ D h A x u z e dx- I XeD h {V x -'VyU 1 )z e dx+ I \\7 X ■ (V(x)\7 x u )z e dx 
Jai Jn" h Jn% u 

- eD h v{- {V x u + V y u 1 )z t ds 

Jv 

- (/ u - (DiVyV )do-)z e dx + e / Div e ■ \7 y voz e ds 

Jn € h JdB(x) JT' 

+ e / XedtU!Z e dx-e D h \7 x m ■ \7 ( Xe z e ) dx 

- e 2 / DiA x v w £ dx - e I D ( (V x ■ \7 y v + \7 y ■ V x v )w e dx + e 2 D t u e ■ V x v w € ds 

+ e 2 / D lV e • (eVw e - V y v ) Ul ds - (eD h A XeUl + 2eD h V Xe ■ V x Ui + 2D h V Xe ■ V y ui)z e dx 

h 

- I dq-q t )-Vu z e dx-e f X eUiq e • Vz e dx, 
Jni v Jnt 



where we have also used that 4> — -0 — z t — w e — eu\ on T e and <j> = eu\ on 
We know that there exist (3 > and 7 > such that 

/3|W|ffi(n«) < / (A l Vz e -Vz e + g £ -Vz e )dx + 7 ||z e || L 2 ( ^ ) , 

J Of 



and we use this to estimate 



iUlWI^fie) + |5 t || We || 2 2(nr) + ^II^IIhi^) + ^Di\\Vw e \\l 2{ni) < j\\z € \\ L 2 {nv +h+I 2 + ... + I n 
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where 
h = 
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h = 



h = 



h = 



h = 



h = 



{ D h W x -{W x u + W v u 1 )z e dx- f \^/ x ■ (V(x)V x u )z € dx 
Jni Jn> h V 

^rtD h J vi ■ (V x u + VyUi)z t ds , 

/ (/ u ■ (DiWyVo) da)z e dx - e / Div Q -V ' y v z e ds 

Jn' h JdB( x ) JT' 

/ (1 - Xe)D h (\/ x ■ V y u 1 )z t dx 
Jni 

\D h (y 2 + 0(e)) • (y x u - V u tti) + D x v\ ■ V y v^jz t ds 
\ (eD h AxeUi + 2eD h \7xe ■ V^ui + 2D h V\t • ^ v ui)z e dx 

e / D h V x U! ■ V(xeZ e ) dx 
Jni 

h 

e 2 / DiA x vow £ dx + e / Di(V x ■ V y Vo + V y ■ V x Vq)w 6 dx — e 2 / Div e ■ V x vow e ds 
Jn\ Jnf Jr* 

e 2 J Div € ■ (eVi> e — \7 y Vo)ui ds 



hi = 



[ i^Q-Qe) -Vu z e dx , 

e / XeUiq e ■ Vz e dx . 
Jni 

h 

For Ii we use that u\ — MVuq and that T>(x) := Dh J Y ( x )(I + ^y^) ^V- We set 

Q(x, y) = V x • ((J + V a M)V x n ) - ^V, • / (I + V y M)V x u dy, 

Jy( x ) 



p(x,y) = -v x ■ (J + V y M), 
and use Lemma I 



5.2 



to obtain I\ < eC||z e ||#i(n|)- Lemma 5.1 
of Lemma |5.2| are satisfied for these choices of Q and p. 
For I 2 we also apply Lemma [5~2"1 this time for the choice 



asserts that the conditions 



1 



V{ x ) JdB( x ) 

p(x,y) = v ■ VyV , 



vq ■ WyVo da, 
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and we again get I 2 < eC\\z e \\ H i^j with C independent of e. For I 3 we have I 3 < 



V^^\\ z A\H 1 (Q' h )- Application of the regularity results in Lemma 



results in I4, + Jg < Ce||£ e ||#i(fje y With the use of Lemma 5.3 and the properties of 



4.8 



and Theorem 



4.9 



X € we estimate J5 + I7 < v / eC||z e ||jyi(f2^)- Another application of the regularity results 
gives I s < eC\\w e \\ H i(n*) and 7 9 < Ce||u ||_f/2 (a ,j(||w e || H 2( O ,) + ||«o||fr»(n;ff2(B(a:))))' For 
Iio and In we use the assumptions on q and q e (as stated in (A2) and (B2)) to get 

ho + hi < y/eC\\z e \\H^si' h y 
Now we obtain 

^d t \\z £ \\ 2 L2(nv + ~$|K||£a ( n r) +D h \\Vz £ \\ 2 L2(nv +e 2 D l \\Vw t \\ 2 L2m < \\z e \\ LHn * h )+ 
eCi||wo||fl- 2 (n^)(ll w e||H 2 (nf) + \\ v o\\H 2 (n-,H 2 (B(x)))) + V^C2\\z e \\ H i^ h - ) + eC 3 \\w e \\ H -L^ 



Using the energy bounds (5.1| and (5.2) together with a Gronwall-type argument lead 
to 

ll z e|lL~(/,L 2 (fi«)) + ll w e|li°=(7,L 2 (^)) + II z e II \ 2 {lM l (U' h )) + 

+ 4we\\h(i,m{ciD) < eC i + V^C 2 \\z e \\ L 2 (ItH i (n%)) , (5.6) 

which, in particular, implies 



\Ze\\ 



L 2 (/,i/!(0' 



■)) < + Vf-C 2 \\z e \\ L 2( I H l( n <-jy 



and thus 



INI^Wo-)) < v%(Cf + VC*2 + 4Ci) 



Combining this with (5.6 1 gives the result 

II 2 £ II L 00 (I,L 2 (tt' h )) + ||™e|U°°(J, L 2 (a«)) + II z e II L 2 (/,if 1 (n«))+e||w e II L 2 (IM 1 < "V^ 

where the constant c is independent of e. The last step uses the evident estimate ||eui(l 
XfOllff 1 ^) ^ C\/ e j and the theorem is proven. 



□ 
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